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We explore heat transfer in molecular junctions between two leads in the absence of a finite net 
thermal bias. The application of an unbiased, time-periodic temperature modulation of the leads 
entails a dynamical breaking of reflection symmetry, such that a directed heat current may emerge 
(ratchet effect). In particular, we consider two cases of adiabatically slow driving, namely (i) periodic 
temperature modulation of only one lead and (ii) temperature modulation of both leads with an ac 
driving that contains a second harmonic, thus generating harmonic mixing. Both scenarios yield 
sizeable directed heat currents which should be detectable with present techniques. Adding a static 
thermal bias, allows one to compute the heat current-thermal load characteristics which includes the 
ratchet effect of negative thermal bias with positive- valued heat flow against the thermal bias, up to 
the thermal stop-load. The ratchet heat flow in turn generates also an electric current. An applied 
electric stop-voltage, yielding effective zero electric current flow, then mimics a solely heat-ratchet- 
induced thermopower ("ratchet Seebeck effect"), although no net thermal bias is acting. Moreover, 
we find that the relative phase between the two harmonics in scenario (ii) enables steering the net 
heat current into a direction of choice. 

PACS numbers: 05.40.-a,44.10.+i,63.22.-m,05.60.Gg 



I. INTRODUCTION 

In recent years, we have witnessed the development 
of nano-devices based on molecular wires [11, [3, H, 0j ■ 
One of their essential features is that the electric cur- 
rent through them can be controlled effectively. One ap- 
proach to such transport control is based on conforma- 
tional changes of the molecule [1, 0, @| • Another scheme 
relies on the dipole interaction between the molecular 
wire and a tailored laser field [1, [H [H [H, [H. A 
further approach employs gate voltages acting on the 
wire [3, UM, [3 ■ The latter allows a transistor- like con- 
trol which has already been demonstrated experimentally 
[T9| . It is therefore interesting to explore as well 
the related control of heat transport. 

In general, heat transport through a molecular junc- 
tion involves the combined effect of electron as well as 
phonon transfer processes. Control of phonon transport 
is much more complicated since the phonon number is 
not conserved. Nevertheless, the field of phononics, i.e. 
control and manipulation of phonons in nanomaterials, 
has emerged [l^l- This includes functional devices, such 
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as thermal diodes [21|, |2^, |23|, |2J, |2 
transistors J2p, [2§| , thermal logic gates 
memories [3l| based on the presence of a static temper- 
ature bias. The corresponding theoretical research has 
been accompanied by experimental efforts on nanosys- 
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tems. In particular, solid-state thermal diodes have been 
realized with asymmetric nanotubes [s^ and with semi- 
conductor quantum dots 

Upon harvesting ideas from the field of Brownian mo- 
tors [s^, m, m, [33, ISll — originally devised for par- 
ticle transport — a classical Brownian heat engine has 
been proposed to rectify and steer heat current in non- 
linear lattice structures [s^ \^ . In the absence of any 
static non-equilibrium bias, a non-vanishing net heat 
flow can be induced by unbiased, temporally alternat- 
ing bath temperatures combined with nonlinear interac- 
tions among neighboring lattice sites. This so obtained 
directed heat current can be readily controlled to re- 
verse direction. If, in addition, a thermal bias across 
the molecule is applied, a heat current can then typi- 
cally be directed even against an external thermal bias. 
This setup is therefore rather distinct from adiabatic and 
nonadiabatic electron heat pumps which involve photon 
assisted transmission and reflection pro cesses in presence 
of irradiating photon sources [1, [ill l42j . 

In this work, we investigate the possibility of steering 
heat through a molecular junction in the presence of a 
gating mechanism. In doing so, the bath temperatures 
of adjacent leads are subjected to slow, time-periodic 
modulations. Both the electronic and the phononic 
heat current are considered, as is sketched in Fig. [T] 
A finite directed ratchet heat current requires break- 
ing refiection symmetry. This can be achieved by spa- 
tial asymmetries in combination with non-equilibrium 
fiel ds 13^. Issl [sgI [stI . [38| or in a purely dynamical way 
[43I |44 l45l kg. |47|| . In this work we focus on an unbiased 
temporal temperature variation in the connecting leads. 

This paper is organized as follows: In Section II, we 
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FIG. 1: Setup of a molecular junction whose electronic level 
El can be gated, while the vibrational frequency u>i is fixed. 
The lead temperatures 2l(r) (t) are subjected to time-periodic 
modulations. 



specify the physical assumptions and introduce our model 
together with the basic theoretical concepts for direct- 
ing heat current across a short, gated molecular junction 
formed by a harmonically oscillating molecule. The heat 
flux is induced by temperature modulations in the con- 
tacting leads. Section III presents the results for case (i) 
where the temperature is modulated in one lead only. We 
elaborate on the phenomenon of pumping heat against a 
static thermal bias and consider the resulting thermo- 
electric power. In Section IV, we consider case (ii) with 
both lead temperatures periodically, but asymmetrically 
modulated. A finite directed heat current emerges from 
harmonic mixing of different frequencies, which entails 
dynamical symmetry breaking. Since the leading order 
of this heat current is of third order only in the driving 
strengths, the overall rectification is weaker as compared 
to case (i) where the heat flux starts out at second or- 
der in the driving amplitude. The latter scheme with its 
symmetric static parameters, however, provides a more 
efficient control scenario: The direction of the resulting 
heat current can be readily reversed either by a gate volt- 
age or by adjusting the relative phase shift within the 
harmonic mixing signal for the temperature modulation. 
Section V contains a summary and an outlook. 



II. PHYSICAL ASSUMPTIONS, MODEL, AND 
BALLISTIC HEAT TRANSFER 

We consider a molecular junction between two leads 
and a static gate voltage acting on the junction. Heat 
transport from both electrons and phonons is taken into 
account. Since we focus on coherent transport in a 
short molecular wire [i^, H^, , electron-phonon inter- 
action can be ignored. Moreover, we treat ballistic heat 
transfer for the electron system and the phonon system 
in the absence of anharmonic interactions and dissipa- 
tive intra-wire scattering processes. Then, the heat flux 
can be obtained in terms of a Landauer-type expression 
involving corresponding temperature-independent trans- 
mission probabilities for both the electrons [Sll, and 
the phonons [13, [E^ . The total Hamiltonian can thus be 
separated into electron and phonon part, i.e.. 



(1) 



each of which consisting of a wire contribution, a lead 
contribution, and an interaction term, such that 
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The short molecular wire is modeled as a single en- 
ergy level and one harmonic phonon mode only. Then 
the Hamiltonian of the wire electron in tight-binding ap- 
proximation reads 



El , 



(3) 



where Ei describes the on-site energy of the tight-binding 
level which can be shifted via a gate voltage. The elec- 
trons in the leads are modeled as ideal electron gases, 
i.e.. 
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where cj^ creates an electron in state \lq) of lead I — L,K. 
The electron tunneling Hamiltonian 



CRqCO+h.C. , (5) 



establishes the contact between the wire and the leads. 
This tunneling coupling is characterized by the spec- 
tral density r/(e) — ^irj^q l^iql'^H^ ~ ^iq)- We assume 
symmetric coupling within a wide-band limit such that 

r,(e) = r. 

The phonon mode is represented by a harmonic oscil- 
lator with the Hamiltonian 
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where Q and P denote the position and the momentum 
operator, respectively, M denotes the atom mass and iOi 
the characteristic phonon frequency of wire. The phonon 
bath and its bilinear coupling to the wire system is de- 
scribed by 
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where xik,pik,^ik are the position operators, momentum 
operators, and frequencies associated with the bath de- 
grees of freedom; mj are the masses and gi = gi, = 
.9r — 9 represent a symmetric phonon wire-lead cou- 
pling strength for lead I — L, R. The position and mo- 
mentum operators can be expressed in terms of the cre- 
ation and annihilation operators for phonons as xik = 
y/h/2miUJik{ali^ + aik) and pik = iyJhmiUJik/2{a\f, - aik)- 
Throughout the following we assume that slow, time 
modulated temperature system acting on the baths 
are always sufficiently slow so that a thermal quasi- 
equilibrium for the molecular wire system can assumed. 
The heat transport then obeys the adiabatic, exact co- 
herent quantum transport laws as discussed in the next 
subsection, see Eq. (IT51) and (HH) below. 
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A. Adiabatic modulation of the lead temperatures 

At thermal equilibrium with temperature T = 
Tl = 7r with equal electro-chemical potentials /il = 
/iR = jjL, the density matrix for the leads read pi cx 
^-{Hf'+Ht-,,iNi)/k^Ti ^ ^Yieve Ni = Y.q4q^iq is the num- 
ber of electrons in lead Z = L, R and k-oTi denotes the 
present lead temperature multipled by the the Boltzmann 
constant. To induce shuttling of heat, we invoke a non- 
equilibrium situation via an adiabatically slow temper- 
ature modulation Ti{t) in the leads. The latter can be 
realized experimentall y, f or example, by use of a heat- 
ing/cooling circulator ^57|. Then the expectation values 
of the electron and phonon lead operators then read 



(alk'aik) = ni{ujk,Ti{t))6wSkk' 



(8) 
(9) 



where Me,Ti{t)) = [exp((e - in) / ^{1)) + 1]^^ and 
ni{uj,Ti{t)) = [exp(?icj/fcBTK^)) - denote the Fermi- 
Dirac distribution and the Bose- Einstein distribution, re- 
spectively, which both inherit a time-dependence from 
the temperature modulation. This implies a time-scale 
separation which is justified by the fact that laser heating 
of a metallic systeni, the electrons undergo rather fast 
thermahzation [sl, [s^ [13, [Ml- The corresponding re- 
laxation times stem from electron-electron and electron- 
phonon interaction, and for a typical metal is in the or- 
der of a few fs or ps, respectively |62i]. Therefore, 
the changes of the lead temperatures occur on a time- 
scale much smaller than the thermal fluctuations itself, 
i.e. 27r/n » 1 ps. 

The time- varying lead temperatures Ti^{t) and Tr(<) 
are assumed to be time-periodic TL(R-)(t) = TL(-R-)(t + 

27r/fi), where Tq = J\Jf) = T^it) denotes the time- 
averaged environmental reference temperature. This im- 
plies a vanishing temperature bias 



Ar(t) = n{t) - Tnit) = 



(10) 



In the long-time limit, the time-dependent, asymptotic 
heat current Jqit) ^ J^{t) + J^^{t) assumes the period- 
icity 27r/r2 of the external driving field 



(11) 



Henceforth we focus on the stationary heat current Jq 
which follows from the average over a full driving period: 
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(12) 



If the lead temperatures are modulated slowly enough 
(adiabatic temperature rocking), the dynamical thermal 
bias AT(i) can be viewed as a static bias at time t in 
the adiabatic limit f2 — > 0. Thus the asymptotic electron 
and phonon heat currents Jq^^^^ (t) can be expressed by 



the Landauer-type formula for electron heat flux [5l|, [S^] 
and for the phonon heat flux [s^, [3, [E^ [E^ , such that 

x[f{e,n{t))~ f{e,Tnm, (13) 
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where T°'(£) and Tp^(ijj) denote the temperature inde- 
pendent transmission coefficients for electrons with en- 
ergy e and phonons with frequency u scattered from left 
lead to right lead, respectively. Note that the two op- 
posite heat fluxes are not at equilibrium with each other 
and that the heat energy transferred by a single electron 
scattering process is e — fi rather than e [4l| . The reason 
for this is the following. At zero temperature, where the 
energy levels below Fermi energy fi are fully occupied, no 
heat current is transferred since no electron can tunnel. 
At finite temperatures, the tunneling process is thermally 
activated. An electron with energy e tunneling from left 
lead to right lead will dissipate to the Fermi energy level. 
Therefore the heat energy transferred by this electron is 
e — fi. 

The electron transmission coefficient T°'(£) can be ex- 
pressed by the electron Green's functions Q 

r'='(e) = Tr[Gt(e)r«G(e)rL], (15) 

where F' = |l)Fi(l|, stems from the tunnel coupling to 
lead Z = L, R. For the present case of a one-site wire, this 
operator is simply a 1 x 1-niatrix, so that the Green's 
function reads 
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(16) 



G(e) = 

wherer = i(rL-frR). 

For a molecular wire with a single level such as de- 
scribed by Eq. ([3]) , the electron transmission assumes the 
Breit-Wigner form and obeys Q 
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(17) 



where we have assumed symmetric electron wire-lead 
coupling such that F = Fl = Fr,. 

The phonon transmission coefficient TP'^(aj) is evalu- 
ated following Ref. |55]. As is shown in Appendix A, 
it assumes for one phonon mode as well a Breit-Wigner 
form, i.e., the temperature-independent phonon trans- 
mission probability reads 



TP'Xt^) 



(18) 



Here wd is the Debye cut 



where 7(w) = ae 
off frequency of phonon reservoirs in the lead and a = 
irg^ / AmMoj'^ incorporating the phonon wire- lead cou- 
pling g = gi^ = g^. 
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B. Experimental parameters and physical time 
scales 

In our numerical investigation we insert the elec- 
tron wire-lead tunnel rate T = 0.11 eV, which has 
been used also to describe electron tunneling between 
a phenyldithiol (PDT) molecule and gold contact [sl ]. 
The phonon frequency uji ~ 1.4 x lO^'^s^^ is typical for 
a carbon-carbon bond 65]. For the Debye cut-off fre- 
quency for phonon reservoirs we use the value for gold 
which is wd = 2.16 x 10^^ s~^. The phonon coupling fre- 
quency a = 1.04 X lO^'"^ s~^ is chosen such that the static 
thermal conductance assumes the value 50 pWK""'^ which 
has been measured in experiments with alkanc molecular 
junctions [66l |. 

These parameters imply physical time scales which 
are worth being discussed. During the dephasing time, 
electron-phonon interactions within the wire destroy the 
electron's quantum mechanical phase. If this time is 
larger than the dwell time, i.e., the time an electron 
spends in the wire, the electron transport is predomi- 
nantly coherent [4^. Following Ref. [i^, we estimate 
the dwell time by the tunneling traversal time r ~ 
h[{Ei — /i)^ -I- r^]~^/^ which for our parameters is of the 
order r 5 fs and, thus, much shorter than the typical 
electron-phonon relaxation time (dephasing time) which 
is of the order of Ips. This implies that the electronic 
motion is predominantly coherent, so that the electron- 
phonon interaction within the wire can be ignored. The 
phonon relaxation time within wire can be estimated as 
1/a 1 fs. Among the above mentioned time scales, the 
maximum time scale is the electron-phonon relaxation 
time which is in the order of ps. Thus the regime of va- 
lidity of our assumption for adiabatic temperature mod- 
ulations is justified when the angular driving frequency is 
much slower the electron-phonon relaxation rate within 
lead, i.e. f7 < 1 THz. 

III. PUMPING HEAT VIA SINGLE-SIDED 
TEMPERATURE ROCKING 

Let us first consider the case in which the temperature 
of one lead is modulated sinusoidally, while the temper- 
ature of the other lead is constant, 

mt) ^To + Acos{nt), (19) 
Tnit) = To. 

Here A and O are the driving amplitude and (angular) 
frequency of the temperature modulation, respectively, 
while To is the reference temperature. The driving am- 
plitude A is positive and bounded by the temperature 
To since Ti^(t) has to remain positive at any time. The 
temperature difference between left and right lead reads 
AT{t) = Acos(r2<), such that the net thermal bias van- 
ishes on time-average, AT{t) — 0. 

The cycle-averaged heat fluxes, both the electronic and 
the phononic one, follow from numerically evaluating the 
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FIG. 2: (Color online) Directed electronic heat current Jq 
as function of onsite energy Ei — for different reference 
temperatures To and temperature oscillation amplitude A = 
30 K. The arrow marks the onsite energy Ei—fi — 0.138 eV for 
which the pumped electron current assumes at temperature 
To — 300 K its maximum. The adiabatic rocking frequency is 
Q. = 3.92 GHz. 
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FIG. 3: (Color online) Total electronic and phononic time- 
averaged, directed heat current Jq as function of the squared 
driving amplitude A'^ with reference temperature at To = 
300 K for the onsite energy Ei - fi = 0.138 eV. The dotted 
line represents the electronic contribution, the dashed line the 
phononic one. The inset depicts the directed heat current as 
a function of the reference temperature To for the amplitude 
A = 30K (yl^ = 900 K^) marked by the arrow in the main 
panel. 



integrals in Eqs. and As expected for an adi- 

abatic theory, we observe that the average heat current 
Jq is independent of the driving frequency f2. This is in 
accordance with the findings for ballistic heat transfer in 
the adiabatic regime. 

In an experiment, the molecular level Ei can be ma- 
nipulated by a gate voltage which influences only the 
electrons. This allows one to tune the electron transport 
while keeping the phonons untouched. In Fig. [2l we de- 
pict the net electron heat current Jq as a function of 
El — fi for a fixed driving amplitude. We find that the 
heat current possesses an extremum for Ei — fi = 0, i.e., 
when the onsite energy is aligned with the Fermi energy. 
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FIG. 4: (Color online) Total directed heat current Jq as the 
function of static thermal bias AT for different driving am- 
plitude strengths A for the temperature modulation. The 
reference temperature is set as To = 300 K and the electronic 
wire level is set as Ei ~ jj, — 0.138 eV. 
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FIG. 5: (Color online) The time-averaged directed electric 
current J<=' as function of the static voltage bias AV for 
different temperature amplitudes A. The reference tem- 
perature is To — 300 K and the electronic wire level is at 
£'i-/i = 0.138eV. 



Interestingly enough, this extremum is a maximum for 
low reference temperature Tq and turns into a minimum 
when the temperature exceeds a certain values. This im- 
plies that the net electron heat current is rather sensitive 
to the on-site energy with respect to the Fermi energy. 
This property thus provides an efficient way to determine 
experimentally the Fermi energy of the wire as an alter- 
native to, e.g., measuring the thermopower as proposed 
in Ref. [6J|. For large gate variations we find that the 
directed electron heat current is significantly suppressed 
since the wire level is far off the electron thermal energy, 
i.e.. El — jjL ^ k^TQ. The directed heat current then 
is dominated by the phonon heat flux. As temperature 
is increased, the peak positions of the pumped electron 
heat current shifts outwards, away from the Fermi energy. 
At room temperature T — 300 K, the peak positions are 
located at Ei ~ = ±0.138eV. 

A. Scaling behavior for small driving strengths 

Figure [3] shows the total heat current Jq as a function 
of the driving amplitude A for the reference temperature 
To = 300 K and the electronic site above the Fermi level. 

For weak driving {A <C Tq), we find Jq''^^^ oc for 
both the electronic and the phononic contribution. This 
behavior can be understood from a Taylor expansion of 
the Fermi-Dirac and the Bose-Einstein distribution, 

5(C,TL)-.g(^,To) 

= 5(C,ro + AT)-5(e,To) ^20) 

= g'{^,To)AT{t) + ^"^Y"^ [AT{t)f + ■ ■ ■ , 

where g represents the Fermi-Dirac function / and the 
Bose-Einstein function n, while g' and g" denote deriva- 
tives with respect to temperature. Note that the time- 
dependence stems solely from the temperature differ- 



ence AT(i). After a cycle average over the driving pe- 
riod, the first term in the expansion vanishes owing to 
AT{t) — 0. Therefore, the leading term of the heat 
current is of second order, i.e. cx [AT(i)]2, which yields 

Jg*-^^-* oc J^^^^[AT{t)]'^dt cx A^, as observed numeri- 
cally. We also plot the directed heat current as the func- 
tion of the reference temperature Tq in the inset of Fig. O 
The directed phonon heat current decreases monoton- 
ically upon increasing the reference temperature. How- 
ever, the emerging total heat current exhibits a relatively 
flat behavior in a large temperature range. This is due 
to the combined effect from phonons and electrons. At 
high temperatures, the electron heat flux dominates the 
overall directed heat flow. 



B. Thermal load characteristics and 
ratchet-induced thermoelectric voltage 

Thus far we have studied heat pumping in the absence 
of a static temperature bias, i.e. for AT{t) = 0. We 
next introduce a static thermal bias such that a thermal 
bias AT := AT{t) ^ emerges. The resulting total 
directed heat current Jq is depicted in Fig. 01 Within 
this load curve, we spot a regime with negative static 
thermal bias AT < and positive-valued overall heat 
flow until AT reaches the stop-bias value, i.e., we find 
a so-called Brownian heat-ratchet effect [3^, |40]. This 
means that heat can be directed against a thermal bias 
from cold to warm like in a conventional heat pump. The 
width of this regime scales with the driving amplitude A^, 
cf. Fig. H 

As can be deduced from Fig. |4l a zero-biased temper- 
ature modulation generates a finite net heat fiow at zero 
temperature bias similar to the heat fiow that would be 
induced by a static thermal bias. 

Near equilibrium, i.e. within the linear response 
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FIG. 6: (Color online) Directed electron heat current Jg as 
function of the wire level Ei— jj. for various reference temper- 
atures Tq. The driving parameters are Ai — A2 = 30K and 
V5 = 0. 



regime, Onsager symmetry relations for transport of con- 
jugated quantities are expected to hold. Therefore, the 
adiabatic temperature modulations are expected to in- 
duce an electric current as well. This net adiabatic elec- 
tric pump current can be obtained by means of the cycle- 
averaged Landauer expression which reads explicitly: 

J'^' = ^ y dt- J de Tie) [fie, n) - /(e, Tr)] . 

(21) 

We in addition apply a net static voltage bias AV. 
Figure \5\ depicts the net electric current- voltage charac- 
teristics J°'(AF) in the presence of an unbiased temper- 
ature modulation while, importantly no external thermal 
bias is applied. For a positive bias voltage AV > 0, the 
net electric current is negative, i.e. the system effectively 
acts as an "electron pump" . 

The value of this externally applied electric stop- 
voltage AVst, which renders the electric current van- 
ishing, mimics here a sole heat-ratchet induced ther- 
mopower. In the present context, this constitutes a novel 
phenomenon which we term ratchet Seebeck effect. Know- 
ingly, the usual thermopower (Seebeck coefficient) is de- 
fined by means of the change in induced voltage per unit 
change in applied temperature bias under conditions of 
zero electric current pJ]. Here, in the absence of a net 
thermal bias, we introduce instead a Griineisen-like re- 
lation, reading 7 — \AVcs/ J'^^\, where AVcs denotes the 
effective static voltage bias which yields the identical elec- 
tric heat current J°' as generated by our imposed temper- 
ature modulation. It is found that this effective voltage 
bias precisely matches the above mentioned stop- voltage, 
i.e. AV^ff — AVst- Due to the linear J'^^-AV character- 
istics, as evidenced with Fig. [SJ the Griineisen-like con- 
stant 7 is independent of the amplitude of the tempera- 
ture modulation. In doing so, we find for Tq = 300 K and 
i?i — /i = 0.138 eV, that this very Griineisen-like constant 
becomes 7 = 61.9 x lO^V/A. 



IV. TEMPERATURE ROCKING IN BOTH 
LEADS: PUMPING HEAT BY DYNAMICAL 
SYMMETRY BREAKING 

We next consider temperature modulations applied to 
both leads in the absence of a thermal bias. The temper- 
ature driving consists of a contribution with frequency f2 
and a second harmonic with frequency 2f2. This entails a 
dynamical symmetry breaking, namely harmonic mixing 
O, S, m, El, 1131 . The time-dependent lead tempera- 
tures are chosen as 

Tl,r = To ± [Ai cos(m) - A2 cos(2m + ip)], (22) 

such that again TYJfj = %Jfj = Tq and AT{t) = 
2[Ai cos{n,t) + A2 cos{2ilt + cp)]. Then the average tem- 
perature bias vanishes irrespective of the phase lag ip. 

In Fig. [HI we depict the resulting electron heat current 
Jq as a function of the on-site energy Ei — ^ for various 
reference temperatures Tq. At low temperatures, the net 
electron heat current exhibits a minimum at the Fermi 
energy. With increasing reference temperature Tq this 
minimum then develops into a local maximum with two 
local minima in its vicinity. The arrow in Fig. [6] marks 
the minimum at Ei - ^ = 0.049 eV for To = 300 K. It is 
interesting that the direction of the net electron heat cur- 
rent can be tuned by the gate variation. For an electron 
wire level close to the Fermi energy, the directed electron 
heat current is negative. By tuning the gate voltage, the 
heat current undergoes a reversal and becomes positive 
when El — fi is larger than 0.15 eV (at reference temper- 
ature To = 300 K) and eventually approaches zero again 
for large detuning. 

Figure [7] shows the corresponding sum of electron and 
phonon heat flow, i.e., the net heat current Jq as a func- 
tion of wire level Ei — ji. The net phonon heat current 

Jq^ is negative for these parameters (not depicted) and is 
not sensitive to the gate voltage. As a consequence, this 
sum of phonon and electron heat transport, Jq, exhibits 
multiple current reversals as the onsite energy Ei — fj, in- 
creases. For small values ofEi—fi, i.e., close to the Fermi 
surface, both the electron and the phonon heat fluxes are 
negative and in phase with the driving. The absolute 
value of the total heat current assumes its maximum (at 
which the heat current is negative). For intermediate 
values of Ei — fi, the direction of total net current Jq 
is reversed due to the dominating positive contribution 
of the electrons. At even larger values of Ei — fi, the 
electron heat current almost vanishes, so that the total 
heat current is dominated by the negative-valued contri- 
bution of the phonons. In the limit of large Ei — /i, we 
find saturation at a negative value. 

We also study with Fig. [5] the net electron and phonon 
heat current as the function of driving amplitudes Ai, 
A2 and the relative phase ip. Both contributions scale 

as Jq^"^^ oc AlA2C08{ip), which implies that they can 
be manipulated simultaneously. This behavior can be 
understood by again expanding the Fermi-Dirac and the 
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FIG. 7; Total net heat current Jq as function of the wire 
level El— jj. for reference temperature To — 300 K, amplitudes 
Ai = A2 = 30 K, and phase lag ip = 0. 

Bose-Einstein function at Tq: 

= To + AT/2) - g{^, To - AT/2) .^S) 
= g\tTo)AT{t) + ^ °' [AT{t)]' + ■ ■ ■ , 

where g = f,n represents the distribution function for 
the electrons and for the phonons, respectively. Note that 
the terms of even order in AT vanish owing to the anti- 
symmetric temperature modulation. Thus, the heat cur- 
rents are governed by the time-average of the odd powers 
[Ar(t)]2"+i with n > 1, since AT = 0. It can be easily 
verified that all these time-averaged odd moments van- 
ish if either amplitude, Ai or A2 vanishes. Note that 
the lowest lowest-order contribution is the third moment 
[AT{t)]^ = 8AlA2Cos{(p). Thus, for small driving am- 
plitudes, Ai, A2 <^ To, the net electron and phonon heat 
current are expected to be proportional to A\A2Costf 
as is corroborated with the numerical results depicted in 
Figs.[51[a,b,c). This proportionality to cos (/?, see Fig.|SlJc), 
is even more robust than a priori expected; this is so be- 
cause the cycle averaged 5-th and 7-th moment are pro- 
portional to If, as well, i.e., [AT(t)]5, [AT{t)]'' oc cos((p). 
This behavior can be employed for a sensitive control 
of the heat current: The direction of the heat current 
can be reversed by merely adjusting the relative phase (f 
between the two harmonics. Note that for the parame- 
ters used in the figure, the net electron heat current Jq 

exceeds the net phonon heat current Jq^ roughly by a 
factor 5. 



V. CONCLUSIONS 

We have demonstrated the possibility of steering heat 
across a gated two-terminal molecular junction, owing 
to lead temperatures that undergo adiabatic, unbiased, 
time-periodic modulations. In a realistic molecule, the 
heat flow is carried by the electrons as well as by the 
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FIG. 8; (Color online) Heat current Jq as function of (a) 
fundamental driving strength Ai , where ^2 = 30 K and (p — 
and (b) as function of the second-harmonic amplitude A2, 
where Ai = 90 K and (fi = 0. (c) Dependence on the relative 
phase (p for Ai = A2 = 30 K. The reference temperature is 
To = 300 K, while the wire level is _Bi - = 0.049 eV. 

phonons. Our study considers both contributions. Two 
scenarios of temperature modulations have been investi- 
gated, namely directed heat flow (i) induced by periodic 
temperature manipulation in one connecting lead and 
(ii) created by a temperature modulation that includes 
a contribution oscillating with twice the fundamental 
frequency. In both cases, we predict a finite heat cur- 
rent which is related to dynamical breaking of reflection 
symmetry. A necessary ingredient is the non-linearity 
of the initial electron and phonon distribution, which is 
manifest in the Fermi-Dirac distribution and the Bose- 
Einstein distribution. The first scenario yields sizable 
heat currents proportional to the squared amplitude of 
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the temperature modulation. The resulting heat flow oc- 
curs in the absence of a static thermal bias. We also stud- 
ied heat pumping against an external static thermal bias 
and computed the corresponding thermal heat-current 
load characteristics. Moreover, the ratchet heat flow in 
turn generates also an electric current. This ratchet heat 
current induces a novel phenomenon, namely a ratchet- 
induced, effective thermopower, see in Fig. [5] 

When the asymmetry is induced by temperature rock- 
ing at both leads, the resulting net heat current becomes 
smaller in size. This is so because the leading-order time- 
averaged heat flow now starts out with the third moment 
of the driving amplitude. The benefit of this second sce- 
nario is the possibility of controlling efficiently both the 
magnitude and the sign of the net heat flow. For example, 
the direction of the heat current can be readily reversed 
via the gate voltage or the relative phase between two 
temperature modulations that are harmonically mixed. 
When adjusting the gate voltage, the directed heat cur- 
rent experiences multiple current reversals. The directed 
heat flow is even up to 7-th order in the amplitude pro- 
portional to the cosine of the phase between the funda- 
mental frequency and the second harmonic. This allows 
robust control even for relatively large temperature am- 
plitudes. 

These theoretical findings may also inspire experimen- 
tal efforts to steer heat in a controlled manner across 
a molecular junction as well as the development of new 
concepts for measuring system parameters via their im- 
pact on the heat current. For example, as elucidated in 
Sect. Ill, the Fermi energy can be sensitively gauged in 
this way. 
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Eq. (33) of that work and assuming symmetric coupling, 
i-e- 7fe,fe'(wo) = if^k'i^o) = 7fe,fc'(^o), Eq. (33) of Ref. ^ 
can be simplified to read 



X! \/--^27fc,fc'(tJo)^fc'('^o) = \l'^Vo,kal, 



k'^k 

(Al) 

where luq is a dummy variable. 

Since we only consider one phonon mode, i.e. k = 1. 
The second term in the left hand side of last equation 
vanishes such that 



[ujI - cJo + 2itJo7i,i(^o)]^i(^o) = \ — Vo,iaJ. (A2) 

V "^0 



Substituting Eq. (46) of Ref. [53], i-e. 



Aki^o) = AkiuJo)Vo,kaL — 
V "^0 



(A3) 



into Eq. ((A2|l . we find 



Ai{ujo) = 



ujf - u!q+ 2ia;o7i,i(wo) 



(A4) 



For one phonon mode, the phonon transmission is de- 
fined from Eq. (48) in 55]. However, this definition is 
l/27r times smaller than the commonly used definition of 
Ref. [13] and (HH]. With the commonly used definition, 
the phonon transmission can be expressed as 



T{uj)^4uj^^l,{uj)\Ai{io)\-\ 



(A5) 



Substituting Eq. (jA4|) into the last equation and omitting 
the subscript in 71,1, we obtain 



APPENDIX A: THE DERIVATION OF PHONON 
TRANSMISSION (fTsjl 



r(c.) 



(A6) 



We derive the phonon transmission coefficient of 
Eq. along the fines of Ref. [SI]. Starting with 



which is the phonon transmission (jlSp employed in the 
main text. 
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